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SIMULTANEOUS DETERMINATION OF FORMATION ANGLES AND 
ANISOTROPIC RESISTIVITY USING MULTI-COMPONENT INDUCTION 

LOGGING DATA 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0001]The invention is related generally to the field of interpretation of measurements 
made by well logging instruments for the purpose of determining the properties of 
earth formations. More specifically, the invention is related to a method for 
simultaneous determination of formation angles and anisotropic fonnation resistivity 
using multi-component resistivity data. 

2. Backgroimd of the Art 

[0002] Electromagnetic mduction and wave propagation logging tools are commonly 
used for determmation of electrical properties of formations surrounding a borehole. 
These logging tools give measurements of apparent resistivity (or conductivity) of die 
formation that when properly interpreted are diagnostic of the petrophysical 
properties of the formation and the fluids therein. 

[0003] The physical principles of electromagnetic induction resistivity well logging 
are described, for example, m, H. G. Doll, Introduction to Induction Logging and 
Application to Logging of Wells Drilled witii Oil Based Mud, Journal of Petroleum 
Technology, vol. 1, p.l48. Society of Petroleum Engineers, Richardson Tex. (1949). 
Many improvements and modifications to electromagnetic mduction resistivity 
instruments have been devised since publication of the Doll reference, supra. 
Examples of such modifications and improvements can be found, for example, in U.S. 
Pat No. 4,837,517; U.S. PaL No. 5,157,605 issued to Chandler et al, and U.S. Pat 
No. 5,452,761 issued to J3e^^ et al. 

[0004] A limitation to die electromagnetic induction resistivity well logging 
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instruments known in the art is that they typically include transmitter coils and 
receiver coils wound so that the magnetic moments of these coils are substantially 
parallel only to the axis of the instrument Eddy currents are induced in the earth 
formations from the magnetic field generated by the transmitter coil, and in the 
induction instruments known m the art these eddy currents tend to flow in ground 
loops which are substantially perpendicular to the axis of the mstrument; Voltages are 
then induced m the receive coils related to the magnitude of the eddy currents. 
Certain earth formations, however, consist of thin layers of electrically conductive 
materials interleaved with thin layers of substantially non-conductive material. The 
response of the typical electromagnetic induction resistivity well logging instrument 
will be largely dependent on the conductivity of the conductive layers when the layers 
are substantially parallel to the flow path of the eddy currents. The substantially non- 
conductive layers will contribute only a small amount to die overall response of &e 
mstrument and therefore their presence will typically be masked by the presence of 
the conductive layers. The non-conductive layers, howev^, are the ones which are 
typically hydrocarbon-bearing and are of the most interest to the instrument user. 
Some earth formations which might be of commercial interest thmfore may be 
overlooked by interpreting a well log made using the electromagnetic induction 
resistivity well logging instruments known m &e art. 

[0005] The effect of formation anisotropy on resistivity logging measurements have 
long been recognized. Kunz and Moron studied the anisotropic effect on the response 
of a conventional logging device in a borehole perpendicular to the bedding plane oft 
thick anisotropic bed. Moron ond Gianzero extended this work to accommodate an 
arbitrary orientation of the borehole to the bedding planes. 

[0006] Rosthal (U.S. Pat No. S329,44S) discloses a method for determining the 
horizontal and vertical conductivities from a propagation or induction well logging 
device. The method assumes that 6, flie angle between the bordiole axis and the 
normal to the bedding plane, is known. Conductivity estimates are obtamed by two 
methods. The first method measures the attenuation of the amplitude of the received 
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signal between two receivers and derives a first estimate of conductivity from this 
attenuation. The second method measures the phase difference between the received 
signals at two receivers and derives a second estimate of conductivity from this phase 
shift. Two estimates are used to give the starting estimate of a conductivity model 
5 and based on this model, an attenuation and a phase shift for the two receivers are 

calculated. An iterative scheme is then used to update the initial conductivity model 
until a good match is obtained between the model output and the actual measured 
attenuation and phase shift. 

10 [0007] Hagiwara shows that the log response of an induction-type logging tool can 
be described by an equation of the form 

F oc ^ (1 - ikL) + iJdie''^ -e^)) 

^ (1) 

where 

JS^ = cos^^+sin^^ 
IS and 

where X is the spacing between the transmitto: and receiver, k is the wavemunber of 
the electromagnetic wave, \i is the magnetic permeability of the medium, 0 is the 
deviation of the borehole angle from the normal to flie formation, X is the anisotropy 
20 &ctor for the formation^ o) is the angular frequency of the electromagnetic v^ave, 0^ is 
the horizontal conductivity of the medium and 6^ is the horizontal dielectric constant 
ofthemediuHL 

[0008] Eq, (3) is actually a pair of equations, one correspondmg to the real part and 
25 one corresponding to the imaginary part of the measured signal, and has two 

unknowns. By makmg two measurements of the measured signal, the parameters k 
and P can be determmed. The two needed measurements can be obtained from (1) R 
and X signals from induction logs, (2) phase and attenuation measurements from 
induction tools, (3) phase or attenuation measurements from induction tools with two 
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different spacings, or (4) resistivity measurements at two different frequencies. In the 
low frequency limit, 6 can be neglected in Eq, (3) and from known values of co and n, 
the conductivity o can be detennmed from k assuming a value of |x equal to the 
permittivity of fiw space 

[0009] Wu (U.S. Patent 6,092,024) recognized that the solution to eqs. (1)- (3) may 
be nonunique and showed how this ambiguity in the solution may be resolved using a 
plurality of measurements obtained from multiple spacings and/or multiple 
frequencies. 

[0010] Strack et al. (U.S. Patent 6,147,496) describe a multicomponent logging tool 
composing a pair of 3-component transmitters and a pair of 3-component receivers. 
Using measurements made at two different frequencies, a combined signal is 
genemted having a reduced dqpendency on the electrical conductivity m Ihe wellbore 
region. 

[0011] United States Patent 5,999,883 issued to Gtpta et al , (the "Ciupta patent"), 
the contents of Msdiich are frilly incorporated here by reference, discloses a method for 
detennination of an initial estimate of the horizontal and vertical conductivity of 
anisotropic earth formations. Electromagnetic induction signals induced by induction 
transmitters oriented along three mutually orthogonal axes are measured. One of the 
mutually orthogonal axes is substantially parallel to a logging instrument axis. The 
electromagnetic induction signals are measured using first receivers each having a 
magnetic moment parallel to one of the orthogonal axes and using second receivers 
each having a magnetic moment perpendicular to a one of the orthogonal axes which 
is also perpendicular to the instrument axis. A relative angle of rotation of the 
perpendicular one of Ihe ordiogoiud axes is calculated from the receiver signals 
measured perpendicular to the instrument axis. An intermediate measurem^t tensor 
is calculated by rotatmg magnitudes of the receiver signals through a negative of &e 
angle of rotation. A relative angle of inclination of one of the orOiogonal axes which 
is parallel to the axis of the mstrument is calculated, from the rotated magnitudes. 
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with respect to a direction of tiie vertical conductivity. The rotated magnitudes are 
rotated through a negative of the angle of inclination. Horizontal conductivity is 
calculated from the magnitudes of the receiver signals after the second step of 
rotation. An anisottopy parameter is calculated from the receiver signal magnitudes 
aftw the second step of rotation. Vertical conductivity is calculated from the 
horizontal conductivity and the anisotropy parameter. One drawback in the teachings 
of Gtipta et al is that the step of determination of the relative angle of inclination of 
the requires measurements of three components of data wifli substantially identical 
transmitter-receiver spacings. Because of limitations on the physical size of the tools, 
this condition is difficult to obtain; consequently the estimates of resistivities are 
susceptible to error. 

[0012] There is a need for a method of simultaneous determination of formation 
angles and anisotropic resistivity that does not require measurements made with 
matched transmitter-receiver spacings. Such a method should prefoabfy 
computational^ effident The present invention satisfies this need. 

SUMMARY OF THE INVENTION 
[0013] The present invention is a metiiod for the simultaneous inversion of 
measurements made by a multicomponent logging tool to obtain a layered resistivity 
model and formation inclination angle and azunuth. A model that includes horizontal 
and vertical resistivities is used to generate a simulated tool response. An iterative 
solution tiiat gives an improved match between the model output and flie field 
observations is obtained using a global objective fimction. The global objective 
function is defined as a sum of a data objective fimction (difference between the 
model oulput and the observed data) and a model otgective function that stabilizes the 
inversion procedure by placing a poialty on large changes in the model at each 
iteration. The measurements may be made by an electromagnetic logging tool havii^ 
an axis mdined to the normal to the bedding planes of the fomiation. The logging 
tool mcludes transmitters and/or recdvei^ with coils mclined to the axis of the tool. 
In a preferred embodim^ of the mvention, the data objective fimction is defined m 
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the coil coordinate system. Surveying of the borehole and orientation sensors on the 
tool provide the necessary information for rotating the model output to the coil 
coordinate system. 

5 BRIEF DESCRIPTION OF T HE FTCTIRTCS 

[00141 The present invention is best understood with refeience to the following 

figures in which like numbers refer to like components. 

FIG. 1 shows an induction instrument disposed in a wellbore penetrating earth 

formations. 

1 0 FIG. 2 shows the arrangement of transmitter and receiver coils in a preferred 
embodiment of the present invention marketed under the name 3DExplorer ™ 
FIG. 3a and 3b shows the angles of relevance in a multicomponent logging tool in an 
inclined borehole. 

FIG. 4 shows the relation between the formation coordinates and the tool and coil 
15 coordinates. 

FIG. 5 is a schematic illustration of the model parameters in an embodiment of the 
inventions. 

FIG. 6 how the position of a straight trajectory in space is described by two angles. 
Fig. 7 shows two sequential rotations of the coils of the logging tool. 
20 Fig. 8 shows the relative dip angle between the formation (Z-O and the trajectory (Z 
Fig. 9 shows the relative rotation angle between the formation {Y^ and the trajectory 

Fig. 10 shows a model used to illustrate the method of the present inventioin and the 
result of inversion of the model output using prior art methods with an incorrect 
25 inclination angle. 

FIG. 11 shows a comparison between the model and the results of using the present 

invention on simulated measurements. 

Fig. 12 shows the residual error as a function of azimuth. 

30 PETAH.EB DESCRIPTION OF THE INVENTION 
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[0015] Referring now to Fig. 1, an electromagnetic induction well logging instrument 
10 is shown disposed in a wellbore 2 drilled through earth formations. The earth 
formations are shown generally at 4. The instrument 10 can be lowered into and 
withdrawn fiom flie wellbore 2 by means of an armored electrical cable 6 or similar 
conveyance known in the art. The mstrument 10 can be assembled from three 
subsections: an auxiliary electronics unit 14 disposed at one end of the instrument 10; 
a coil mandrel unit 8 attached to the auxiliary electronics unit 14; and a 
receiver/signal processing/telemetry electronics unit 12 attached to the other end of 
the coil mandrel unit 8, this unit 12 typically being attached to the cable 6. 

.i^[0016] The coil mandrel unit 8 includes induction transmitter and receive coils, as 
will be fiirther explained, for inducmg electromagnetic fields in the earth formations 4 
and for receiving volta^ signals mduced by eddy currents flowing in the earth 
formations 4 as a result of the electromagnetic fields induced tfierein. 

[0017] The auxiliary electronics unit 14 can include a signal generator and power 
amplifiers (not shown) to cause alternating currents of selected fiiequencies to flow 
through transmitter coils in the coil mandrel unit 8. 

[0018] The receiver/signal processing/telemetry electronics unit 12 can include 
receiver circuits (not shown) for detecting voltages induced in receiver coils in the 
coil mandrel unit 8, and circuits for processing these received voltages (not shown) 
into signals representative of the conductivities of various layers, shown as 4A 
flirough4Fofthe earth formations 4. As a matter of convenience flie receiver/signal 
processing^l«netry electronics unh 12 can indude signal telemetry to transmit the 
conductivity- related signals to the earth's sut&ce along the cable 6 for further 
processing, or alternatively can store the conductivity related signals in an ^tpropriate 
recording device (not shown) fat processmg after the instrument 10 is withdrawn 
fix>m the wellbore 2. 

[00191 Referring to Fig. 2, the configuration of transmitter and receiver coils in a 
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preferred embodiment of the 3DExplorer ™ induction logging instrument of Baker 
Hughes is shown. Three orthogonal transmitters 101, 103 and 105 that are referred to 
as the Tx, T^, and Ty transmitters are shown (tixe z- axis is tibie longitudinal axis of the 
tool). Corresponding to the transmitters 101, 103 and 105 are associated receivers 
107, 109 and 111, referred to as the R,, R^, and Ry receivers, for measuring the 
corresponding magnetic fields. In a preferred mode of operation of the tool, the 
Hyy and H^^ components are measured, though other components may also 

be used. 

[0020] There are four coordmate systems involved in the processing of 
multicomponent induction logging data, which are the earth, formation, tool, and coil 
coordinate systems, respectively. The earth coordinate system comprises the z- axis 
in the direction of gravity, the x- axis pointing to the magnetic North, and the y- axis 
in the direction of magnetic East The other three coordinate systems are defined in 
the earth coordinate system. Fig. 3a shows the formation coordinate system. Thez- 
axis of the formation coordinate system, Z^, is normal to the bedding plane, vdiile 
and F^are on the bedding plane. The two formation angles, ^and (2^, are to be 
recovered fix>m our inversion processing. The tool coordinate system is described by 
{X\Y\Z^m Fig. 3b. The z-axis of the tool system is normal to flie tool plane and X' 
and 7^ are on the tool plane. The dip and azimuth angles for the borehole trajectory 
are denoted by JFand /'respectively. The data are measured in the coil coordinate 
system QC ^ F Z 0 that shares the same z-axis with the tool coordinate system. 

[0021] The forward modeling method used in this paper is the one developed by Epov 
et al. (1997), which requires the input of relative dip and rotation angles under the 
formation coordinate system (Fig. 4). This requires determination of the relative dip 
and azimuth angles using 9, $,r and 7 as outlined in Appendix A. Oncetiiedata 
are generated under &e formation coordinate system, two rotations axe p^ormed to 
convert the data from formation system to the coil syst^. The magnetic field data are 
first rotated finom the formation coordinate system into an intermediate magnetic 
field,/?', on the tool plane 
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[0022] Any magnetic field component H^ia the formation coordinate system can be 
rotated to an intennediate magnetic field H' by the rotation 



(4) 



where the rotation matrix Rg is given by 



cos^ 0 sin^ 

0 1 0 
- sin^ 0 cos^ 



(5) 



The intermediate magnetic field is fiirther related to the measured magnetic field If 
by 



10 H" = R^H% 
where 



.(6) 



15 



cos^ sin^ 0 
- sin^ cos^ 0 
0 0 1 



(7) 



The fi)rmation angles may then be calculated as (see Appendix B) 



tan2^ = 



and 



HL-HL 



(8) 



2h: 



rrc rrc 



(9) 
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[0023] This is the method used in Gupta and, as noted above, requires that all field 
components involved must be measured at the same spacing. This is difScult to 
achieve in practice, and, in the present invention, the angles (9 and 0 are determined 
simultaneously along with the resistivity parameters. 

5 

[00241 After conversion, the 3D coordinate system is reduced to a 2D system 
described by relative dip 9 and relative aamuth angle ^ .This is shown in in Fig. 5. 
The horizontal axis 151 is the x- axis in the formation coordinate system and the 
vertical axis 153 is the z- axis Z^in the formation coordinate system. The model is a 
10 1-D model with layers 1, 2, ....1, 1+1, the individual layers bemg characterized by 
horizontal and vertical resistivitiesy?^ ^ and a The data D in an induction loggmg 

e}q)eriment can be e3q)ressed as a nonlinear fimction of the model parametm 
nonlinear function 

^=/(w) (10) 
15 where m = (0, 0,p^pj where /> *s are vectors of horizontal and vertical 

resistivities in the model In a preferred embodiment of tiie model, fhe angles 0 and 
0 are defined with respect to an absolute coordinate system where fhe z- axis is fhe 
vertical, and the x- axis is the true north, though any other refi^ence coordinate 
system may be used. 

20 

[0025] The ultimate goal of inv^on is to recover ^and 0 in additional to the 
resistivity structure. There are two possible ways to recover formation angle 
information. The first one is to include the absolute formation angles, 0md , in the 
model objective fimction and recover them directly fiom the inversion. The second 
25 method is to solve for the relative formation dip angle first, and then recover absolute 
formation dip and azimuth angles via a post inversion procedure. 

[002d] The model parameters ni^ be determined by first defimng a d^ 
function : 
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0^ = lWiD'"^-D)f (11) 

where D^^ and D are the observed data and the data predicted from the model 
respectively, and is a suitable weighting matrix. 

5 [0027] Those versed in the art would recognize that attempting an optimization 

procedure with respect an objective function such as that given by eq. (1 1) relative to 
the model parameters may involve numerical instabilities. Prior art methods have 
attempted to reduce these mstabilities by selecting step fectors that are used to update 
the model parameters. The present invention uses a different procedure for this 
10 stabilization by including, in the objective function, quantities that are related to the 
modelitself. It is desirable fliat the model = p/j atthen-th 

iteration be close to flie model w '^^ = (0"-^ 0"-', p^'^ , p^^) at the (n-l )- th 
iteration. 

15 [0028] In one embodiment of the invention, this is done by defining separate 

objective functions are defined for the horizontal and vertical resistivities, formation 
dip and formation azimuth as: 




and 
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In 



n-\ 



(15) 



The use of logaritbmic fimctions ensures non-negative solution of the recovered 
model and allows tiie model parameters for resistivity, dip angles and thicknesses to 
span the same numerical range in the inversion. 



10 



The model objective function is then defmed as 

0„ = a^O^ + a^O, + OgOg + a^O^ (16) 

where the five coefficients ar^ , «^ and ar^are determined from the following 
conditions: 



15 



20 



a,,+a^+ag+a^=l (17) 

i.e., the contributions to ibe cost fiinction fix)m the four components are substantially 
the same at each iteration. 
The solution to eq. (18) is 



with 



(18) 



I = O^O.Og + 0^0,0^ + 0,0^0^ + O.OgO^ (19) 
and similar expressions for the remaining a coiefficients. 

[0029] The four weighting coefficients are calculated and tq>dated at each iteration. 
Tlien for tiie preferred embodiment, 



12 



BNSOCXJID: <WO_0aa7323eA?_L> 



wo 02/073236 



PCT/US02/06832 



5 




and the four weighting coefiBcients are of the form given by eqs. (18) and (19). 
[0030] This gives a model objective function of the form 

10 O^^a^Of^ + a^O^^^a^O^+a^O^ (24) 

A global objective function is then defined as a combination of the data objective 
function firom eq. (10) and the model objective function from eq. (25) as 

15 0=0„-^fi-\0,-0'^) (25) 

vdiere/?*^ is the LagrangianmiUtiplier and is tiie target misfit level (a noise 
level). 
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[0031] The ribnlinear problem of optimizing eq. (26) is done by using linearization i 
an iterative process. The global objective function at the n+i iteration is given by 

0=lw„[Sm^m^''^ -n^f^fi-'ljlV,(D^ ^ f[m^"^]^ JSmf ^ O^j (26) 

where J =[J^Jy,J0,J4^] is the global sensitivity matrix. The sensitivities for 
horizontal and vertical resistivity are generated using an adjoint Green's function 
solution while the sensitivities for the formation angles are calculated using a finite 
difiference metiiod. 

[0032] The global model weighting matrix is given by 





. 0 


0 


0 " 


0 




0 


0 


0 


0 




0 


0 


0 


0 





(27) 



Those versed in the art would recognize tiiat if aU the a's are set equal to 1, 
matrices W,, and are set equal to unitary matrices and if a'*' = 0, then the mversion 
method of the present invention is equivalent to tiie prior art Marquardt-Levenberg 
method For this reason, the'metfaod of the present equation described above may be 
called a generalized Marquardt-Levenberg method. With increasing iterations, the 
global objective function is weighted more and more towards the data misfit temi so 
that when the solution is close to the optimal, the global objective function is 
substantially equal to the data objective function. 

[0033] Even though it is possible to invert the 3DEK data to simultaneously recover 
formation dip and azimuth angles as well as horizontal and vertical resistivity, 
reduction in the number of inversion parameteis is still welcome for reducing the 
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computational time. We may exclude foraiation azimuth angle fix)m the inversion 
scheme using the following relations from Appendix B: 

HI, + = H'„ + 

Hl=K (28) 
Eq. (28) indicates that the sum of the xx- and>y- components of the measured 
magnetic field inside a layered earth is irrelevant to the determination of (p. Further 
more, we note that if^ at ^^AS"". 



[0034] In an alternate embodiment of the mvention, we use HP and (jff /f ^) tl as 
the data in the inversion and fix the rotation angle ^ at 45° , and exclude formation 
azimuth angle as mversion parameter. This is compromised with the exclusion of the 
two cross components from the inversion. Given the fact that the cross components 
Tisually have much lower signal-to-noise ratio than the three main components, their 
1 5 exclusion fi-om ID inversions do not significantly affect the final results. In an 

alternate embodiment of the invention, in addition to the three main components, 
measurements are also made the xy- and yx-components, and use is made of the 
following relationships to include cross components m the inversion (Appendix B): 

20 ^r;-/r; = fr;-^; (29) 

Note that the coil separations for yz- and xy-components must be the same in order to 
make use of eq. (28). The signal level of cross components over ID earfli are usually 
much smaller than those of three main components, and subtraction of these cross 
component m Eq. (29) may fiirth^ reduce the signal to noise ratio of die data. This 
25 means that in order to take advantages of Eq. (29), the data quality for cross 
components must be high. 

[0035] The objective fimction for this alternate embodiment of the invention 
becomes 



15 



BNSOOCIO: <WQ_0207323eA2_L> 



wo 02/073236 



PCT/US02/06832 



0,„ = a;,0^+a^O,+aeO, 



(30) 



where 




(31) 



2 



wdth 



2 = O,O,+O,0e+0eO,(32). 
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After 0 is obtained by the inversion procedure, a post inversion procedure described 
in Appendix C is used to extract information about the absolute formation dip and 
azunuth. 

[0036] In yet another embodiment of the invention, the objective function described 
inamediately above by eqs. (30) - (32) is used wilh flie inversion carried out only for 
^(Fig. 8), Rf^, and A possible range of #is obtained fcom geologic and 
geophysical information and #is digitized into N disorete values 0'^[0i,02» ^3 
.... <P„ ] (Fig. 9). For each of these discrete values of 0, , an mversion is performed 
for 0, , and R^ and the data misfit determined for each of these discrete values. On 
the basis of these mdividual inversions, the model with the least data misfit (smallest 
value of the data objective function) is selected. Those versed in the art would 
recognize that inclusion of the model objective function in the iterative process is 
primarily for the purpose of stabilizing the inversion. 

EXAMPLE 

[0037] The results of using the meflxod of the present invention are illustrated by an 
example. Fig. 10 shows an example ofa resistivity model in which there is a 
constant value of /Rf, of 2.0 as mdicated by 208 in the right track. The left track 
shows the actual values ofR^ 203 and Rf, 205. 
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[0038] Synthetic data were generated for an inclination angle of 60^ The curves 207 
and 205 in Fig. 10 show the results of a mverting the synthetic model output wherein 
an incorrect inclination angle of 55® was assumed and the inversion was limited to 
determination of the resistivity values only, i.e., the inclination angle was not included 
in the mversion process. The inversion was carried out using prior art methods. The 
curve 209 shows the detemiined value of R^/R^ using the incorrect inclination angle. 
It is clear that not including the inclination angle in the inversion process leads to 
large errors in the inversion if a wrong value of the inclination angle is assumed. 

[0039] Fig. 11 shows the result of inverting the synthetic model data generated using 
the model shown in F^* 10 using the present method, i.e,, a simultaneous inversion 
ofresistivity and angles. The initial values for bofli and are 5 ohm-m whole- 
space, and the starting value for relative dip is 75^ The standard deviation is set to 
1% of the data strength plus a negligible constant base value of 10"^ A/m. This 
constant threshold can effectively handle zero-crossings on Hxx, Hyy , Hxy , and Hxz 
components. The inversion converged to the desired misfit level after 9 iterations. 
The recovered resistivity mcKiel in Fig. 11 is a good representation of the true model. 
The anisotropy ratio has been perfectiy resolved. The recovered relative dip is 60.06''. 

[0040] To illustrate how to recover information about relative azimuth angle, we 
calculate tiie x^-misfit as the function of relative azimuth in Fife. 12. The x^-misfit 401 
reaches a min i mum at tiie true solution <H)^ Because the data are contaminated 
with noise, the inversion can only recover a narrow zone around the true solution for 
the relative azimuth angle. The higher the noise level 403, the wider the zone. The 
X^-misfit reaches anoflier minimum at <(»+7C=l 80**, since the magnetic field is 
symmetric over a ID earth. Hence a priori information is needed to determine the 
final solution for the relative azimuth angle. This example indicates that our 
simultaneous inversion algorithm can efiectively handle zero-crossings in the data 
and provide reliable information about the relative angle as well as the anisotropic 
resistivity structures. 
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[0042] The present invention has been discussed above with respect to measurements 
made by a transverse induction logging tool conveyed on a wireline. This is not 
intended to be a lunitation and the method is equally ^plicable to measurements 
made using a comparable tool conveyed on a measurement-while-drilling (MWD) 
5 assembly or on coiled tubing. 

[0043J While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the art. It is 
intended that all variations within the scope and spirit of the appended clainois be 
1 0 embraced by the foregoing disclosure. 

APPENDKA 
Calculate 6. 7} and g in the earth coordinate RY stftm 

15 [0042] The borehole trajectory is measured in the field, and at any depth alone Ae 
borehole the azunufli T of flie borehole and its inclmation and T are known. The 
dkectional vector for the Z-axis of tiie tool coordinate system is 

Z' = sinr cos^r + sinr sin^J + cosFit (A-1) 
where 

20 I J, and ^are the directional unit vectors along Ae X-, Y-, and Z^lnections m the 
earth coordinate system. Similariy, we have for the formation coordinate system 

= sin® cos of + sin© sin^ + cos®^ (A-2) 

25 The relative dip between and 2? is 

^= arccos(Z-^ • Z') (A-3) 

Because tiie fonnation and tool coordinate systems share the same Y-axis, i.e., Y^Y*, 
therefore 
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= r =: X r (a-4) 

where x denotes the vector cross product 

5 [0043] The relative azimuth angle is calculated using the measured relative bearing, 
which is the angle between and high side vector of the tool plane. The high side 
vector, h, is the vector in the tool plane pointing to the high side of the tool. Vector h 
can be computed in the following manner: 

10 [0044] When r is smaller than or equal to 90^ 

h = sin(r + 90^)cosOF)/ + sin(r + 90^)smOF)7 + cos(r + 90°)jt (A-S) 



and Mdien r is greater than 90*^ , 

h = sin(270^ - F) cos(^ + 1 80^ )i + sin(270*^ - F) sinC^' + 1 80^ ) J + cos(270° - r)jt 
15 (A-6) 
The angle betwera h and is then 

K= arccos(r' - A) (A-?) 

The relative angle is given by 

<*=^+90°-v. ^ 

20 These two angles, 9zsA ^ are then used in the forward modeling to generate 
magnetic field in the fonnation coordinate system. 



APPENDIX B 

25 The Exclus ion of Relative Rotation angle firom mversion 

[0045] For forward modeling purposes, we need to convert the magnetic field from 
the formation coordinate system to the coil coordinate system. The magnetic field 
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in the fonnation coordinate system can be rotated to an intermediate magnetic field 

H': 



H' = R^H^R^ 

where the rotation matrix is given by 



(B-1) 



COS0 0 sin^ 

0 1 0 
- sin^ 0 cos^ 



(B-2) 



and 



ffi K Hi 



xz 



K K K 
hL K Hi 



(B-3) 



M„ tl^ 



K 

HL 



(B-4) 



1 0 This intermediate magnetic field is fiirther related to tiie measm«d magnetic field IF 



H' = RlH% 
■wiiere 



(B-5) 







sin(i( 


0 




-sin^ 




0 




0 


0 


1 



(B-6) 
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The definitions of ^and ^ are given in Fig. 3 and 4. Expanding eq. (B-5) gives 

-sin<i>cos<^' +sin^dH* 
•9' ^ yy 

^f, =cosjJsin(*/r' -sin^^' 
xy ^ ^ XX ^ yx 

xy ^ yy 

xz jcy yz 

yx XX yx 



xy ^ yy 



(B-7) 



yy ^ XX ^ ^ yx 

yz xy yz 

7x ^ zx zy 

2Z 2Z (B-g) 



S [0046] If we closely examine Equations (B-7) and (B-8), we find that 

" " (B-9) 
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Eq. (B-9) is the same as Eq. (29). The amplitudes of cross components are much 
smaller than those of the three main components over a layered earth. Neglectmg 
cross components m Equation (B-7) leads to sunplified angle estimation in Eqns. (8) 
and (9): 



tan2^= 



2flL 



tan2^= 



(B-10) 



10 



[0047] The relative azimuth, can be calculated from the measured data using 
equation (B-IO). We tiien rotate the measured magnetic field into the intermediate 
field If, and conq)ute ^fi:om equation (B-10). 



APPENDIXC 

Determination of Fomuttion Azimuth Angle 



15 



[00481 After the relative formation dip angle is determined, the formation azimuth 
angle can be calculated via post inversion process. To simplify mafli, let us consider a 
straight well trajectory. Mathematically, the three rotations discussed above and 
shown in Figs. 6 and 7 m^ be ejqjressed as follows: 
cos(r^) sin(r„) o" 













L 



0 0 1 

cosT 0 - sinr 

0 1 0 
sinF 0 cosF 



cos^ sinT 0 
-sinT cos^ 0 
0 0 1 



(C-1) 
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where /, , and /. are the unit vectors in the Earth coordinate system; i,' , iy ' and f. ' 
are unit vectors in the tool coordinate system.; and is the tool rotation angle. 



10 



15 



[0049] The possible directions of the formation vector constitute a cone around the 
trajectory (Figf. 9). The angle between the cone gaierator and the trajectory equals ft 
To determine the trae orientation of the formation, the tool rotation angle, Trb, must 
be used in combination with a cross-component measurement The method of 
transforming the inversion results into true dip «& azimuth angles includes the 
following steps: 



1. 
2. 



Select the cone generator m the gtavity/trsgectory plane, above the tr^ectory. 
Position Ae tool sensors as follows: z- along the trajectory; x - in the gravity 
trajectory plane, perpendicular to flie trajectory and upwards. This is called 
the "principal position." 

Calculate the magnetic tensor, defined in tije principal tool system. Note 
that in the principal tool position, all ooss-components are zero. 



20 



XX 

0 



H, 



ZX 



H, 



yy 
0 



XZ 

0 



ZX _ 



(C-2) 



4. Rotate the tensor ^about the trajectory <p* to satisfy the individual 
measurements 

jOrC rjC , rjc 

•» and >y in addition to the abeady satisfied sum ^«"^^xy. Rotation 
of the tensor J7 is described by the foilowmg equation: 
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cosgi sm.f 0 
-sin^ cosj> 0 
0 0 1 



H, 



XX 

0 



zx 



H 



yy 
0 



xz 

0 



a JL 



cos^ -sin^ 0 
smp cosf 0 
0 0 1 



(C-3) 



[00501 We can rotate H by the angle T^^ to bring the sensors jc and >; from the 
principal position to the actual tool position. After the rotation, the sum of new 
components 

5 H^-hH'yy 

will have a correct value (as the sum is independent of the rotation angle) but the 
individual components wiU not have correspond to the actual measwements. To 
satisfy the measurements, the formation vector must be rotated about the trajectoiy. 

10 [0051] Instead of rotating the formation vector, one can rotate the tool in the opposite 
directioiL Given this feet, the following algorithm may be used to satisfy each of the 
measured components 



15 



(a) Rotate the tensor by a necessary angle 9 to satisfy the following equation: 



_ cos^ 9H„ +sin^ VH^ jH^ +H^)HHxx -Hy^)cos2p HL 
H'^ sm^ (pH„ +cos2 ipH^ iH^^H^)-(H^ -H^)cos2p " //^ ^ " 

4) 

(b) The solution to eq, (C-4) is not unique. Two angles <p and (p +7C satisfy the 
equation. To resolve the non-uniqueness, the cross component /f^ is used. 
20 From eq. (C-3), tiiis is 
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H;^ = COS ^ (C-5) 



[0052] Comparing the signs of calculated and measured components, from eq.(C-5). 
the correct solution of eq.(C-4) may be established. It is worth noting that the spacing 
5 for the component differs from that of components H„ and H^. Theoretically, to 
obtain a correct synthetic value for H'je,eq. (C-3) must be recalculated with the 
correct spacing of the . NevCTtheless, if we use only the sign of the H'^ 
component, eq. (C-3) may be sufficient with a single spacing. 

.10 [0053] Another non-uniqueness is related to the direction of Hie tool rotation. Both 
angles, <p and -(p, satisfy Equation (C-4). This non-uniqueness may be resolved using 
cross-component (or HyJ. It follows from Equation (C-3): 



— H + H 

15 =("/^« +^vK)cos(^)sin(?^) = ^sm(2p) 

2 (c-6) 

[0054] The sign of the cross-component allows us to determine the direction of 

the tool rotation. Remarks regarding the spacing of the cross-component H^^ remain 

true for the cross-component H^. Let be the solution to Equations (C-4) and (C-5). 

20 (c) Select the cone generator at the angle (-4>o) and rotate both tool and formation 
about trajectory by the angle T^. Sudi a rotation puts sensors m their real 
positions with respect to the Earth coordinate syst^ At the same time, eq. 
(C-4) remains satisfied because tibe solution depends only on the relative 
position of the formation and s^or. After rotation, the selected cone 

25 generator moves from the angular position (-Oq) to the position <I>=-<I>o+Trb- 

This is the formation vector F. The following equation describes the vector F: 



F = sin(F„;)[cos(0)/; -sm(0)i;)]+cos(F„,)/; 
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Here, vectors i',, and i'^ represent the orientation of the tool sensors in the 
principle positioa They may be calculated fiom Equation (C-1) assuming Frb=0: 



0 1 0 
,sin(r^) 0 cos(rrf,p)J 



-sin(r^) cos(rrf„) 0 

0 0 i; 



y 



5. 



(C-8) 



Project vector F into the gravity system using Equations (C-7) and (C-8): 



= [sin(i^«/)cos(0) -sin(F^)sin(«D) cosC^-^)] 



(C-9) 



6. Calculate the true formation dip F^^ and azimuth fiom eq. (C-9) 



tan(F„) = -^ = 



sm(F^)[cos(<I»)cos(r^)co8(2'^)+sin(<I>)sin(r^)] + cos(F^)sin(r,^)cos(r^) 



10 



and 



cos(F^) = -sin(F;^)cos(<l>)sin(r^)+cos(F^,)cos(rrf^) 



(C-10) 
(C-U) 
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CLAIMS 

WHAT IS CLAIMED IS: 

1 . A method of logging of subsurface fonnations including a plurality of layers 
each having a horizontal resistivity and a vertical resistivity, the method 
comprising: 

(a) conveying an electromagnetic logging tool into a borehole in the 
subsurface formations; 

(b) obtaining multi-component measurements indicative of said horizontal 
and vertical resistivities using said logging tool; 

(c) defining a model of said plurality of layers, said model including a 
horizontal resistivity and a vertical resistivity associated with each of 
said plurality of layers, an azunudi angle and a inclination angle of 
said layers referenced to an earth coordinate system; 

(d) determmmg expected responses of the loggmg tool to said model; 

(e) definmg a data objective function related to a difference between said 
e^qpected responses and said multi-component measurements; and 

(f) using an iterative procedure including stabilization for updating the 
model to reduce said data objective function. 

2. The method of claim 1 wherein the stabilization comprises: 

(a) defining a model objective function related to changes of component 
values of said model in successive iterations; 

(b) defining a global objective function as a weighted sum of the data 
objective function and the model objective function; and 

(c) updating the model to using a sensitivity of tiie global objective 
function to model parameters. 

3. The method of claim 1 wherein said electromagnetic tool has an axis inclined 
to a normal to said layers. 
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1 4. The method of claim 1 wherein obtaining said multi<omponent 

2 measurements further comprises using at least one transmitter coil having an 

3 axis inclined to an axis of the logging tool. 

1 5. The method of claim 1 wherein obtammg the multi-component measurements 

2 further comprises using at least one receiver coil having an axis inclined to an 

3 axis of the loggmg tool. 

1 6. The method of claim 1 wherein obtainmg the multi-component measurements 

2 further comprises using at least two receiver coils having an axis inclined to 

3 an axis of the logging tool. 

1 7. The method of claim 1 wherein said multi-component measurements comprise 

2 Hjo^Hj^and i^L measurements. 

1 8. The method of claim 7 wherein the multicomponent measurements further 

2 comprise and /i^measuremenb. 

1 9. The method of claim 2 wherem the model objective function fiirttier 

2 comprises a weighted sum of mdividual objective functions related to said 

.3 horizontal and vertical resistivities of the plurality of layers, the dip angle and 

4 the azimuth angle. 

1 10. The method of claim 1 vsdierein the model objective function is related to a 

2 ratio of value of at least one component of the model after an iteration to a 

3 value of the at least one component before said iteration. 

1 11. The method of clahn 1 )^ereindetennining said expected responses furth 

2 conqnises: 

3 (i) obtaining an azunuth and mclinationofsaid borehole at a depth; 

4 (ii) obtaining a toolface orientation of said logging tool; 
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5 (iii) detenninmg from said obtained borehole incliM^ 

6 said obtained toolface orientation, a relative azimuth and inclination of 

7 the tool to the formations. 

1 12. A method of logging of subsurface formations including a plurality of layers 

2 each having a horizontal resistivity and a vertical resistivity, the method 

3 comprising: 

4 (a) conveying an electromagnetic logging tool into a borehole in the 

5 subsurface formations; 

6 (b) obtaining multi-component measurements indicative of said horizontal 

7 and vertical resistivities using said logging tool; 

8 (c) defining a model ofsaid plurality oflayers, said model including a 

9 horizontal resistivity and a vertical resistivity associated with each of 

10 said plurality oflayers, and a inclination angle ofsaid layers 

11 referenced to an earth coordinate sfystem; 

^12 (d) determinmg expected responses oftfae logging tool to said model using 

13 a relative azimuth angle between the tool and the formation; 

14 (e) defining a data objective function related to a difference between said 

1 5 e3q>ected responses and said multi-component measurements; and 

16 (f) using an iterative procedure including stabilization for updating the 

1 7 model to reduce said data objective function. 

1 13. The method of claim 12 further comprising summing and components 

2 ofsaid measurements and usmg said summed component in the data objective 

3 function. 

1 14. The mediod of claim 1 2 herein said electromagnetic tool has an axis inclined 

2 to a normal to said layers. 
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The method of claim 12 wherein obtaining said multi-component 
measurements further comprises using at least one transmitter coil having 
axis inclined to an axis of the logging tool. 



1 1 6. The method of claim 12 wherem obtaining the multi-component 

2 measurements further comprises using at least one receiver coil having an axis 

3 inclined to an axis ofthe logging tool. 

1 1 7. A method of logging of subsurface formations including a plurality of layers 

2 each having a horizontal resistivity and a vertical resistivity, the method 

3 comprising: 

4 (a) conveying an electroipagnetic logging tool into a borehole in the 

5 subsurface formations; 

6 (b) obtaining multi-component measurements indicative of said horizontal 

7 and vertical resistivities using said logging tool; 

8 (c) defining a model of said plurality of layers, said model including a 

9 horizontal resistivity and a vertical resistivity associated with each of 

10 said plurality oflayers, an azimuth angle and a inclination angle of 

1 1 said layers referenced to an earth cooidinate system; and 

12 (d) defining a plurality of azimuth angles and, for each of said plurality of 

13 azimuth angles, 

14 (i) mverting said multi-component measurements to give a model 

15 and determining a data misfit for said model, 

16 (ii) deteraiining one of said plurality of azimuth angles for which 

17 said data misfit attains a minimum value. 
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SIMULTANEOUS DETERMINATION OF FORMATION ANGLES AND 
ANISOTROPIC RESISTIVITY USING MULTI-COMPONENT INDUCTION 

LOGGING DATA 

BACKGROUND OF TH E INVENTION 

1. Field of the Invention 

[0001]The invention is related generally to the field of interpretation of measurements 
made by well logging mstruments for the purpose of detennining the properties of 
earth formations. More specifically, the invention is related to a method for 
sunultaneous determination of formation angles and anisotropic formation resistivity 
using multi-component resistivity data. 

2. Background of the Art 

[0002] Electromagnetic induction and wave propagation logging tools are commonly 
used for determination of electrical properties of formations surrounding a borehole. 
These logging tools give measurements of apparent resistivity (or conductivity) of the 
formation that when properly mterpreted are diagnostic of the petrophysical 
properties of the formation and the fluids therein. 

[00031 The physical principles of electromagnetic induction resistivity well logging 
are described, for example, m, H. G. Doll, Introduction to Induction Logging and 
Application to Logging of Wells Drilled with Oil Based Mud, Journal of Petroleum 
Technology, vol. 1, p.l48. Society of Petroleum Engmeers, Richardson Tex. (1949). 
Many improvements and modifications to electromagnetic induction resistivity 
instalments have been devised since publication of the Doll reference, supra. 
Examples of such modifications and improvements can be found, for example, in U.S. 
Pat. No. 4,837,517; U.S. Pat. No, 5,157,605 issued to Chandler et al, and U.S. Pat. 
No. 5,452,761 issued to Beard ct al. 

10004] A limitation to the electromagnetic induction resistivity well logging 
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instruments known in the art is that they typically include transmitter coils and 
receiver coils wound so that the magnetic moments of these coils are substantially 
parallel only to the axis of the instrument Eddy currents are induced in the earth 
formations fix)m the magnetic field generated by the transmitter coil, and in the 
induction instruments known in the art these eddy currents tend to flow in ground 
loops which are substantially perpendicular to the axis of the instrument. Voltages are 
then induced in the receiver coils related to the magnitude of the eddy currents. 
Certain earth formations, however, consist of thin layers of electrically conductive 
materials interleaved with thin layers of substantially non-conductive material. The 
response of the typical electromagnetic induction resistivity well logging instrument 
will be largely dependent on the conductivity of the conductive layers when the layers 
are substantially parallel to the flow path of Ae eddy currents. The substantially non- 
conductive layers will contribute only a small amount to the overall response of the 
instrument and therefore their presence will typically be masked by the presence of 
the conductive layers. The non-conductive layers, however, are the ones which are 
typically hydrocarbon-bearing and are of the most interest to the instrument user. 
Some earth formations vAnch might be of commercial interest therefore may be 
overlooked by interpreting a well log made using the electromagnetic induction 
resistivity well logging instruments known in the art 

[0005] The effect of formation anisotropy on resistivity logging measurements have 
long been recognized. Kunz and Moran studied the anisotropic effect on the response 
of a conventional logging device in a borehole perpendicular to the bedding plane oft 
thick anisotropic bed. Moran and Gianzero extended this work to accommodate an 
arbitrary orientation of the borehole to the beddmg planes. 

[0006] Rosthal (U.S. Pat No. 5,329,448) discloses a method for determining the 
horizontal and vertical conductivities fiom a propagation or induction well logging 
device. The method assumes tfiat 6, the angle between the borehole axis and the 
normal to the bedding plane, is known. Conductivity estimates are obtained by two 
methods. The first method measures the attenuation of the amplitude of the received 
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signal between two receivers and derives a first estimate of conductivity from this 
attenuation. The second method measures the phase diflFerence between the received 
signals at two receivers and derives a second estimate of conductivity from this phase 
shift. Two estimates are used to give the starting estimate of a conductivity model 
5 and based on this model, an attenuation and a phase shift for the two receivers are 
calculated. An iterative scheme is then used to update the initial conductivity model 
until a good match is obtamed between the model output and the actual measured 
attenuation and phase shift. 

10 [0007] Hagiwara shows that the log response of an mduction-type loggmg tool can 
be described by an equation of the form 



Voijj{- 2e^^ (1 - ikL) + ikl{e^^ - e"^ )) 



(1) 



where 

>2 ^ 2 



(2) 



yff' = cos' ^+ sin' ^> 
IS and 

where L is the spacing between the transmitter and receiver, k is the wavenumber of 
the electromagnetic wave, |i is the magnetic permeability of the medium, 6 is the 
deviation of the borehole angle fix>m the normal to the formation, X is the anisotropy 
20 factor for the formation, o) is the angular firequency of the electromagnetic wave, o^ is 
the horizontal conductivity of the medium and is the horizontal dielectric constant 
of the mediimi. 

[0008] Eq. (3) is actually a pair of equations, one corresponding to the real part and 
25 one corresponding to the imaginary part of the measured signal, and has two 

unknowns. By making two measurements of the measured signal, the parameters k 
and p can be determined. The two needed measurements can be obtained from (1) R 
and X signals from induction logs, (2) phase and attenuation measurements from 
induction tools, (3) phase or attenuation measxirements from induction tools with two 
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different spacings, or (4) resistivity measurements at two. different frequencies. In the 
low frequency limit, € can be neglected in Eq. (3) and from known values of (o and |i, 
the conductivity a can be determined from k, assuming a value of \i equal to the 
permittivity of free space 

5 

[0009] IVu (U.S. Patent 6,092,024) recognized that the solution to cqs, (1)- (3) may 
be nonunique and showed how this ambiguity in the solution may be resolved using a 
plurality of measurements obtained from multiple spacings and/or multiple 
frequencies. 

10 

lOOlOJ Stracketal. (U.S. Patent 6,147,496) describe a multicomponent logging tool 
comprising a pair of 3-component transmitters and a pair of 3-component receivers. 
Using measurements made at two different frequencies, a combined signal is 
generated having a reduced dependency on the electrical conductivity in the wellbore 
IS region. 

[00111 United States Patent 5,999,883 issued to Gupta et al , (the "Gupta patent"), 
the contents of which are fully incorporated here by reference, discloses a method for 
determination of an initial estimate of the horizontal and vertical conductivity of 

20 anisotropic earth formations. Electromagnetic induction signals induced by induction 
transmitters oriented along three mutually orthogonal axes are measured. One of the 
mutually orthogonal axes is substantially parallel to a logging instrument axis. The 
electromagnetic induction signals are measured using first receivers each having a 
magnetic moment parallel to one of the orthogonal axes and using second receivers 

25 each having a magnetic moment perpendicular to a one of the orthogonal axes which 
is also perpendicular to the instrument axis. A relative angle of rotation of the 
perpendicular one of the orthogonal axes is calculated from the receiver signals 
measured perpendicular to the instrument axis. An intermediate measurement tensor 
is calculated by rotating magnitudes of the receiver signals through a negative of the 

30 angle of rotation. A relative angle of inclination of one of the orthogonal axes which 
is parallel to the axis of the instrument is calculated, from the rotated magnitudes, 

4 
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with respect to a direction of the vertical conductivity. The rotated magnitudes are 
rotated through a negative of the angle of inclination. Horizontal conductivity is 
calculated from the magnitudes of the receiver signals after the second step of 
rotation. An anisotropy parameter is calculated from the receiver signal magnitudes 
after the second step of rotation. Vertical conductivity is calculated from the 
horizontal conductivity and the anisotropy parameter. One drawback in the teachings 
of Gupta et al is that the step of determination of the relative angle of inclination of 
the requires measurements of three components of data with substantially identical 
transmitter-receiver spacings. Because of limitations on the physical size of the tools, 
this condition is difScult to obtain; consequently the estimates of resistivities are 
susceptible to error. 

[0012] There is a need for a method of simultaneous determination of formation 
angles and anisotropic resistivity that does not require measurements made with 
matched transmitter-receiver spacings. Such a method should preferably 
computationally efficient. The present invention satisfies this need. 

SUMMARY OF THE INVENTION 
[0013] The present invention is a method for the simultaneous inversion of 
measurements made by a multicomponent logging tool to obtain a layered resistivity 
model and formation inclination angle and azimuth. A model that includes horizontal 
and vertical resistivities is used to generate a simulated tool response. An iterative 
solution that gives an improved match between the model output and the field 
observations is obtained using a global objective ftmctioa The global objective 
Amotion is defined as a sum of a data objective function (di£ference between the 
model output and the observed data) and a model objective function that stabilizes the 
inversion procedure by placing a penalty on large changes in the model at each 
iteration. The measurements may be made by an electromagnetic logging tool having 
an axis inclined to the normal to the bedding planes of the formation. The logging 
tool includes transmitters and/or receivei? with coils inclined to the axis of the tool. 
In a preferred embodiment of the invention, the data objective function is defined in 
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the coil coordinate system. Surveying of the borehole and orientation sensors on the 
tool provide the necessary information for rotating the model output to the coil 
coordinate system. 

BRIEF DESCRIPTION OF TH E FIGURES 
[0014] The present mvention is best understood with reference to the following 
figures in which like numbers refer to like components. 
FIG. 1 shows an induction instrument disposed in a wellbore penetrating earth 
formations. 

FIG. 2 shows the arrangement of transmitter and receiver coils in a preferred 
embodiment of the present invention marketed under the name 3DExpiorer ™ 
FIG. 3a and 3b shows the angles of relevance in a multicomponent logging tool in an 
mclined borehole. 

FIG. 4 shows the relation between the formation coordinates and the tool and coil 
coordinates. 

FIG. 5 is a schematic illustration of the model parameters in an embodiment of the 
inventions. 

FIG. 6 how the position of a straight trajectory in space is described by two angles. 
Fig. 7 shows two sequential rotations of the coils of the logging tool. 
Fig. 8 shows the relative dip angle between the formation (Z^ and the trajectory (Z 0- 
Fig. 9 shows the relative rotation angle between the formation (K^ and the trajectory 

(n. 

Fig. 10 shows a model used to illustrate the method of the present inventioin and the 
result of inversion of the model ou^ut using prior art methods with an incorrect 
inclination angle. 

FIG. 11 shows a comparison between the model and the results of usmg the present 

invention on simulated measurements. 

Fig. 12 shows the residual error as a function of azimuth. 

DETAILED DESCRIPTION OF THE INVENTION 
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[0015] Referring now to Fig. 1, an electromagnetic induction well logging instrument 
10 is shown disposed in a weilbore 2 drilled through earth formations. The earth 
formations are shown generally at 4. The instrument 10 can be lowered into and 
withdrawn from the weilbore 2 by means of an annored electrical cable 6 or similar 
5 conveyance known in the art. The instrument 10 can be assembled from three 

subsections: an auxiliary electronics unit 14 disposed at one end of the instrument 10; 
a coil mandrel unit 8 attached to the auxiliary electronics unit 14; and a 
receiver/signal processing/telemetry electronics unit 12 attached to the other end of 
the coil mandrel unit 8, this imit 12 typically being attached to the cable 6. 

10 

[0016] The coil mandrel unit 8 includes induction transmitter and receiver coils, as 
will be further explained, for inducing electromagnetic fields in the earth formations 4 
and for receiving voltage signals induced by eddy currents flowing in the earth 
formations 4 as a result of the electromagnetic fields induced therein. 

15 

[0017] The auxiliary electronics unit 14 can include a signal generator and power 
amplifiers (not shown) to cause alternating currents of selected frequencies to flow 
through transmitter coils in the coil mandrel unit 8. 

20 [0018] The receiver/signal processing/telemetry electronics unit 12 can include 

receiver circuits (not shown) for detecting voltages induced in receiver coils in the 
coil mandrel unit 8, and cux:uits for processing these received voltages (not shown) 
into signals representative of the conductivities of various layers, shown as 4A 
through 4F of the earth formations 4. As a matter of convenience the receiver/signal 

25 processing/telemetiy electronics unit 12 can include signal telemetry to transmit the 
conductivity- related signals to the earth's surface along the cable 6 for fiirther 
processing, or alternatively can store the conductivity related signals in an appropriate 
recording device (not shown) for processmg after the instrument 10 is withdrawn 
from the weilbore 2. 

30 

[0019] Referring to Fig. 2, the configuration of transmitter and receiver coils in a 

7 
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preferred embodiment of the SDExpIorer ™ induction logging instrument of Baker 
Hughes is shown. Three orthogonal transmitters 101, 103 and 105 that are referred to 
as the T„ T„ and T, transmitters are shown (the z- axis is the longitudinal axis of the 
tool). Corresponding to the transmitters 101, 103 and 105 are associated receivers 
5 107. 109 and 111, referred to as the R,. and R, receivere, for measuring the 

corresponding magnetic fields. In a preferred mode of operation of the tool, the H„ 
Hyy and H„ components are measured, though other components may also 

be used. 

1 0 [0020] There are four coordinate systems involved in the processing of 

multicomponent induction logging data, which are the earth, formation, tool, and coil 
coordinate systems, respectively. The earth coordinate system comprises tiie z- axis 
in tiie direction of gravity, the x- axis pointing to the magnetic North, and the y- axis 
in the direction of magnetic East. The other three coordinate systems are defined in 

15 the earth coordinate system. Fig. 3a shows the formation cooidinate system. Thez- 
axis of the formation coordinate system, Z^, is normal to tiie bedding plane, while 
and r^are on tiie bedding plane. The two formation angles, ©and are to be 
recovered from our inversion processing. The tool coordinate system is described by 
iX',Y'.Z 0 m Fig. 3b. The z-axis of the tool system is normal to the tool plane and X' 

20 and K' are on the tool plane. The dip and azimutii angles for the borehole trajectoiy 
arc denoted by Fand /^respectively. The data are measured m tiie coil coordinate 
system (A", ¥',2 that Glares the same z-axis witii the tool coordinate syston. 

10021] The forward modeling metirad used in this paper is the one developed by Epov 
25 et al. (1 997), which requires tfie input of relative dip and rotation angles under tiie 

formation coordinate system (Fig. 4). This requires determination of tiie relative dip 
and azunutii angles using 6, F and 7 as outiined in Appoidix A. Once flie data 
are generated under the formation coordinate system, two rotations are performed to 
convert the data from formation system to the coil system. The magnetic field data are 
30 first rotated Scorn the formation coordinate system into an intermediate magnetic 
field, /f', on the tool plane 
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I0022J Any magnetic field component H^m the fonnation coordinate system can be 
rotated to an intomediate magnetic field ^' by the rotation 

H' = R^H^R, (4) 
where the rotation matrix Rgis given by 



Re- 



COS0 0 sin^ 

0 1 0 
-sin^ 0 cos^ 



(5) 



The intermediate magnetic field is fijrther related to the measured magnetic field If 
by 



10 H' = RjH'R^ 
w^ere 



(6) 



15 



cos^ sin^ 0 
-sin^ cos^ 0 
0 0 1 



(7) 



The fonnation angles may then be calculated as (see Appendix B) 



tan 2^ = 



and 



2HL 



(8) 



2H' 

tan2^ = — ^ 



(9) 
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[0023] This is the method used in Gupta and, as noted above, requires that all field 
components involved must be measured at the same spacing. This is difficult to 
achieve in practice, and, m the present invention, the angles &and <Pare determined 
simultaneously along with the resistivity parameters. 

[00241 After conversion, the 3D coordinate system is reduced to a 2D system 
described by relative dip ^and relative azimuth angle ^ .This is shown in in Fig. 5. 
The horizontal axis 151 is the jc- axis A'^in the fomiation coordinate system and the 
vertical axis 153 is the z« axis Z^in the formation coordinate system. The model is a 
1-D model with layers 1, 2, ...J, 1+1, the individual layers being characterized by 
horizontal and vertical resistivities p^^ and p^j. The data Diaan induction logging 
experiment can be expressed as a nonlinear function of the model parameters m via a 
nonlinear function 

D=f(m) (10) 
where m = (&, 0, p^ p^) where p *s are vectors of horizontal and vertical 
resistivities in the model. In a preferred embodiment of the model, the angles ^and 
0 are defined with respect to an absolute coordinate system where the z- axis is the 
vertical, and the x- axis is the true north, though any other reference coordinate 
system may be used. 

[0025] The ultimate goal of inversion is to recover ^and 0 in additional to the 
resistivity structure. There are two possible ways to recover formation angle 
information. The first one is to include the absolute formation angles, &and <P , in the 
model objective function and recover them directly from the inversion. The second 
method is to solve for the relative formation dip angle first, and then recover absolute 
formation dip and azimuth angles via a post inversion procedure. 

[0026] The model parameta:s may be determined by first defining a data objective 
function : 
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^ \\1V(D^'' - D)f (11) 

where and D ait the observed data and the data predicted from the model 
respectively, and is a suitable weighting matrix. 

5 [0027] Those versed in the art would recognize that attempting an optimization 

procedure with respect an objective function such as that given by eq. (1 1) relative to 
the model parameters may involve numerical instabilities. Prior art methods have 
attempted to reduce these instabilities by selecting step factors that are used to update 
the model parameters. The present invention uses a different procedure for this 
10 stabilization by including, in the objective function, quantities that are related to the 
model itself It is desirable that the model m " = (^0", Pa", atthe«-th 
iteration be close to the model w '^^ = <p""', p^'^' , pj^^) at the (n-J)- th 
iteration. 

IS [0028] In one embodiment of the invention, this is done by defining separate 

objective functions are defined for the horizontal and vertical resistivities, formation 
dip and formation azimuth as: 
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The use of logarithmic functions ensures non-negative solution of the recovered 
mode! and allows the model parameters for resistivity, dip angles and thicknesses to 
span the same numerical range in the inversion. 

5 

The model objective function is then defined as 

On, = a,0^ + a,0, + a^O^ + a^O^ (16) 

vAiere the five coeflBcients , and er^are determined fix)m the following 
conditions: 

10 a,0^ = a,0, = a^O, = a^O^ 

ay^ + flf, =1 (17) 

i.e., the contributions to the cost function from the four components are substantially 
the same at each iteration. 
1 5 The solution to eq. (1 8) is 

^ (18) 

with 

Z = O^O^Og + 0^0,0^ + O.OgO^ + O^OfO^ (19) 
and similar expressions for the remaining a coefficients. 

20 

[0029] The four weighting coefficients are calculated and updated at each iteration. 
Then for the preferred embodiment. 



In 



dv 



(15) 



12 
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and the four weighting coefficients are of the fonn given by eqs. (1 8) and (1 9). 
[0030] This gives a model objective function of the fonn 

10 O^^a^Oh+a^O.^a^O^-^a^O^ (24) 

A global objective function is then defined as a combination of the data objective 
function fix)m eq. (10) and the model objective function fiom eq. (25) as 

15 0=0„^fi-\0,-0'^) (25) 

where ^ is the Lagtangian multiplier and (7" is the target misfit level (a noise 
level). 

13 
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(0031) The nonlinear problem of optimizing eq. (26) is done by using linearization in 
an iterative process. The global objective function at the «+7 iteration is given by 

0 = ||fr„[j;n + m^"> ^/i^^l^yff-^l^^ (26) 

where y = [J^ Je^J^] is the global sensitivity matrix. The sensitivities for 
horizontal and vertical resistivity are generated using an adjoint Green's function 
solution while the sensitivities for the formation angles are calculated using a finite 
difiference mediod. 

[0032] The global model weighting matrix is given by 





. 0 


0 


0 


0 




0 


0 


0 


0 




0 


0 


0 


0 


V«7. 



(27) 



Those versed in the art would recognize that if all the a*s are set equal to 1, the 
matrices FF^ and are set equal to unitaiy matrices and if - 0, then the inversion 
method of the present invention is equivalent to the prior art Marquardt-Levenberg 
method. For this reason, die method of the present equation described above may be 
called a generalized Marquardt-Levenberg method. With increasing iterations, the 
global objective function is weighted more and more towards the data misfit temi so 
that when the solution is close to the optimal, the global objective function is 
substantially equal to the data objective function. 

[0033] Even though it is possible to invert the 3DEX data to simultaneously recover 
formation dip and azimuth angles as well as horizontal and vertical resistivity, 
reduction in the number of inversion parameters is still welcome for reducing the 

14 
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computational time. We may exclude fonnation azimuth angle fiom the inversion 
scheme usmg the following relations from Appendix B: 

^==K (28) 
Eq. (28) indicates that the sum of the xx- and yy- components of the measured 
magnetic field inside a layered earth is irrelevant to the determination of 0. Further 
more, we note that ff„= at <J)=45''. 



i0034J In an alternate embodiment of the invention, we use H'„ and /2 as 

the data in the inversion and fix the rotation angle ^ at 45" , and exclude formation 
azimuth angle as inversion parameter. This is compromised with the exclusion of the 
two cross components from the inversion. Given the feet that the cross components 
usually have much lower signal-to-noise ratio than the three main components, their 
exclusion fiom ID inversions do not significantly affect the final results. In an 
alternate embodiment of the invention, in addition to the three main components, 
measurements are also made the xy- and yx-components, and use is made of the 
following relationships to include cross components in the mversion (Appendix B): 

^;-/^; = /f;-/^; (29) 

Note that the coil separations for yz- and jQ'-components must be the same in order to 
make use of eq. (28). The signal level of cross components over ID earth are usually 
much smaller than those of three main components, and subtraction of these cross 
component in Eq. (29) may fiirflier reduce the signal to noise ratio of the data. This 
means that in order to take advantages of Eq. (29), the data quality for cross 
components must be high. 

[0035] The objective fimction for this alternate embodiment of the invention 
becomes 
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0„ = a^O^ + a,0, + aeOe (30) 
where 



with 

5 S = 0,0,+0,Oe+OeO,(32). 



After 8 is obtained by the inversion procedure, a post inversion procedure described 
in Appendix C is used to extract information about the absolute formation dip and 
azimuth. 

10 

[0036] In yet another embodiment of the invention, the objective function described 
immediately above by eqs. (30) - (32) is used with the mversion carried out only for 
^(Fig. 8)» Ri,, and R^. A possible range of <Z>is obtained fiom geologic and 
geophysical information and <2> is digitized into N discrete values 0-[0i,02* ^3. 

IS "'^n] (Pig* I^OT each of these discrete values of (Z^^ , an mversion is performed 
for 0, , and R^ and the data misfit determined for each of these discrete values. On 
the basis of these individual inversions, the model with the least data misfit (smallest 
value of the data objective function) is selected. Those versed in the art would 
recognize that inclusion of the model objective function in the iterative process is 

20 primarily for the purpose of stabilizing the inversion. 

EXAMPLE 



[0037] The results of using flie method of the present invention are illustrated by an 
25 example. Fig. 10 shows an example of a resistivity model in which there is a 

constant value of R^/R^ of 2.0 as indicated by 208 in the right track. The left track 
shows the actual values of R^ 203 and R^ 205. 



16 
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(0038J Synthetic data were generated for an inclination angle of 60". The curves 207 
and 205 in Fig. 10 show the results of a inverting the synthetic model output wherein 
an incorrect inclination angle of 55* was assumed and the inversion was limited to 
determination of the resistivity values only, i.e., the inclination angle was not included 
in the inversion process. The inversion was carried out using prior art methods. The 
curve 209 shows the determined value of R,/R, using the incorrect inclination angle. 
It is clear that not including the inclination angle in the inversion process leads to 
large errors in the inversion if a wrong value of the inclination angle is assumed. 

10039] Fig. 11 shows the result of inverting the synthetic model data generated using 
the model shown in Fig. 10 usmg the present method, i.e., a simultaneous inversion 
of resistivity and angles. The initial values for both p^ and are 5 ohm-m whole- 
space, and the starting value for relative dip is 75*». The standard deviation is set to 
1% of the data strength plus a negligible constant base value of 10 ' A/m. This 
constant threshold can effectively handle zero-aossings on Hxx, Hyy, Hxy, and Hxz 
components. The inversion converged to the desired misfit level after 9 iterations. 
The recovered resistivity model in Fig. 11 is a good representation of the true model. 
The anisotropy ratio has been perfectly resolved. The recovered relative dip is 60.06". 

[0040] To illustrate how to recover information about relative azimuth angle, we 
calculate the x^-misfit as the function of relative azimuth in Fig. 12. The x^-misfit 401 
readies a minimum at the true solution Because the data are contaminated 
wifli noise, the mversion can only recover a narrow zone around the true solution for 
the relative azimuth angle. The higher the noise level 403, the wider the zone. The 
X^-misfit reaches another minimum at ()>+7i=180*', since the magnetic field is 
symmetric over a ID eardi. Hence a priori information is needed to determine the 
final solution for the relative azimuth angle. This example indicates that our 
simultaneous inversion algorithm can effectively handle zero-crossings in the data 
and provide reliable information about the relative angle as well as the anisotropic 
resistivity structures. 
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[00421 The present invention has been discussed above with respect to measurements 
made by a transverse induction logging tool conveyed on a wireline. This is not 
intended to be a limitation and the meAod is equally applicable to measurements 
made using a comparable tool conveyed on a measurement-while-drilling (MWD) 
5 assembly or on coiled tubing. 

[00431 While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the art. It is 
intended that all variations within the scope and spirit of the appended claims be 
1 0 embraced by the foregoing disclosure. 

APPENDIX A 
Calculate 6. 6> Tl and 7} in the earth coordinate system 

1 S [00421 The borehole trajectory is measured in the field, and at any depth alone the 
borehole flie azimuth T of the borehole and its inclination and F are known. The 
directional vector for the Z*axis of the tool coordinate system is 

Z' = sinF cos^r + sinF sin^J + cos Fife (A-1) 
where 

20 i Jy and ^are the directional unit vectors along the and Z-directions in the 
earth coordinate system. Similarly, we have for the formation coordinate system 

Z-^ = sin© cosOr + sin© smOj + cos©^ (A-2) 

25 The relative dip between and Z* is 

^=arccos(Z^ -Z') (A-3) 

Because the formation and tool coordinate systems share the same Y-axis, i.e., ¥^=7*, 
therefore 

18 
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= r = X r (a^) 

where x denotes the vector cross product. 

5 [0043] The relative azimuth angle is calculated using the measured relative bearing, C, 
which is the angle between and high side vector of the tool plane. The high side 
vector, /i, is the vector in the tool plane pointing to the high side of the tool. Vector h 
can be computed in the following manner: 

1 0 [0044] When F is smaller than or equal to 90° 

h = sin(r + 90^)cos(^)/ + sin(r + 90^)sin(^)7 + cos(r + 90^)jt (A-5) 

and when F is greater than 90** , 

/i = sin(270° -r)cos(^ + 180%r + sin(270^ -r)sin(T + 180^)7 + cos(270^ ^F)* 
15 (A-6) 
The angle between h and is then 

V = arccos(7' • A) (a-7) 

The relative angle is given by 

jlS = C+90--V. (A.8) 

20 These two angles, OzxA ^, are then used in the forward modeling to generate 
magnetic field in the formation coordinate system. 



APPENDIX B 

25 The Exclusion of Relative Rotation angle from inversion 

[0045] For forward modeling puiposes, we need to convert the magnetic field from 
the formation coordinate system to the coil coordinate system. The magnetic field W 
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in the formation coordinate system can be rotated to an intermediate magnetic field 



where the rotation matrix is given by 



(B-1) 



COS0 0 sin^ 

0 1 0 
-sin^ 0 cos^ 



(B-2) 



and 



Hi K 

K K K 

Hi K Hi 



(B-3) 



n„ H„ 
^'y ^'yy ^'yt 



Hi 



(B-4) 



} 0 This intermediate magnetic field is further related to the measured magnetic field If 
by 



vhexe 



(B-5) 





cos^ 


sin^ 


0 




-sin^ 


cos^ 


0 




0 


0 


1 



(B-6) 
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The definitions of ^and ^ are given in Fig. 3 and 4. Expanding eq. (B-5) gives 

Hi = cos 2 ^h' - cos ^ sin ^H^ 
XX XX yx 

-sin^ cos + sin ^ <f ^ 
//f, =cos^Jsin^' -sia^dH^ 

xy XX ^ yx 

+ cos ^ + cos ^isin SH' 
xy yy 

xz jcy >!Z 

^f^ = cos^sin<tt?^ +cos^«*^' 
yx ^ ^ XX ^ yx 



and 



-sin^^i,, -cos^sindiff' 
xy yy 



H^yy = sin^ + cos^Jsin#^ 

+sin^cos^' +coP'dH^ . 

xy yy 

yz ^ xy ^ yz 

zx zx ^ zy 

zz zz 



(B-7) 



(B-8) 



[0046] If we closely examine Equations (B-7) and (B-8), we find that 



(B-9) 
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Eq. (B-9) is the same as Eq. (29). The amplitudes of cross components are much 
smaller than those of the three main components over a layered earth. Neglecting 
cross components in Equation (B-7) leads to simplified angle estimation in Eqns. (8) 
and (9): 

2fr' 



tan2^=- 



2H' 

tan2^= ^ 



" ^ (B-IO) 



10 



[0047] The relative azimuth, ^, can be calculated from the measured data using 
equation (B-10). We then rotate the measured magnetic field into the intermediate 
field If, and compute Ofcom equation (B-IO). 



APPENDKC 

Determination of Formation Azimuth Angle 



IS [0048] After the relative formation dip angle is determined, the formation azimuth 

angle can be calculated via post inversion process. To simplify math, let us consider a 
straight well trajectory. Mathematically, the three rotations discussed above and 
shown in Figs. 6 and 7 may be expressed as follows: 

cos(7;j sin(r^) 0 
-sin(r;^)cos(7;jO 
0 0 1 



y 
I 



cosr 0 -sinr 

0 1 0 
sinr 0 cosF 



cosV sin^ 0 
-sinT cos^ 0 
0 0 1 



JL'rJ 



(C-1) 
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where , iy and /. are the unit vectors in the Earth coordmate system; i,' , iy' and /. ' 
are unit vectors in the tool coordinate system.; and Tp^ is the tool rotation angle. 

[0049] The possible directions of the formation vector constitute a cone around the 
5 trajectory (Fig. 9). The angle between the cone generator and the trajectory equals ft 
To determine the true orientation of the formation, the tool rotation angle, Trb, must 
be used in combination with a cross-component measurement. The method of 
transforming the inversion results into true dip & azimuth angles includes the 
following steps: 

10 

1 . Select the cone generator in the gmvity/trajectory plane, above the trajectory. 

2. Position the tool sensors as follows: z- along the trajectory; x - m the gravity 
trajectory plane, perpendicular to the trajectory and upwards. This is called 
the '^principal position.** 

15 3. Calculate the magnetic tensor,//, defined in the principal tool system. Note 
that in the principal tool position, all cross-components are zero. 
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0 



yy 
0 



0 



(C-2) 



4. Rotate the tensor H about the trajectory (p® to satisfy the individual 
measurements 



fjC ffC TJC . wjC 

^ and yy in addition to the ahready satisfied sum yy . Rotation 



of the tensor //is described by the following equation: 
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cosg> sinp 0 
-sin^ cos^ 0 
0 0 1 



H 



XX 

0 



zx 



0 



yy 
0 



cos^ -sin^ 
sin^ cos^ 
0 0 



0 
0 
1 



(C.3) 



[0050] We can rotate ^by the angle Tgg to bring the sensors xwxiy from the 
principal position to the actual tool position. After the rotation, the sum of new 
components 

5 H'„^H'„ 

will have a correct value (as the sum is indq)endent of the rotation angle) but the 
individual components will not have correspond to the actual measurements. To : 
satisfy the measurements, the formation vector must be rotated about the tiajectoiy . 

10 lOOSlJ Instead of rotating the formation vector, one can rotate the tool in the opposite 
direction. Given this feet, the following algorithm may be used to satisfy each of the 
measured componoits 



15 



(a) Rotate the tensor by a necessaiy angle <p to satisfy the following equation: 



H'„ cos^ 9H„ 4sin^ ^Hy, iH„+Hj^)+iH„ -iy^)cos2^ ff^ 
H'y, sin^ <pH„ +cos2 fH„ ' (H„+H^)-(H„ -H^)cos2f> ' fl^ 

4) 

(b) The solution to eq. (C-4) is not unique. Two angles <p and (p +u satisfy the 
equation. To resolve the non-uniqueness, the cross component is used. 
20 From eq. (C-3), this is 
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H'^ =//^COS^(C-5) 

[0052] Comparing the signs of calculated and measured components, from eq.(C-5). 
the correct solution of eq.(C-4) may be established. It is worth noting that the spacing 
for the component H„ differs from that of components H„ and Theoretically, to 
obtain a correct synthetic value for H'^ eq. (C-3) must be recalculated with the 
correct spacing of the H„ . Nevertheless, if we use only the sign of the H'^ 
component, eq. (C-3) may be sufiBcient with a su^le spacing. 

[0053] Anofter non-uniqueness is related to the direction of the tool rotation. Both 
angles, <p and -<p, satisfy Equation (C-4). This non-uniqueness may be resolved using 
cross-component (or HyJ. It follows from Equation (C-3): 



=(-^xr +^;y)cos(?»)sin(^) = ^ — ^s\n{2f) 

^ (C-6) 

[0054] The sign of the cross-component allows us to determine the direction of 
the tool rotation. Remarks regarding the spacing of the cross-component remain 
true for the cross-component H^. Let Oo be Ae solution to Equations (C-4) arid (C-5). 

(c) Select the cone generator at the angle (-Oo) and rotate both tool and formation 
about trajectory by the angle T^. Such a rotation puts sensors in their real 
positions with respect to the Earth coordinate system. At the same time, eq, 
(C-4) remains satisfied because the solution depends only on the relative 
position of the formation and sensor. After rotation, the selected cone 
generator moves from the angular position (-Oq) to the position 0=-Oo+Trb. 
This is the formation vector F. The following equation describes the vector F: 

F = sin(F,,/ )[cos(0)/; - sin(O)/; )] -h cos(iVe/ )i'z 
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Here, vectors i\, i\, and i'^ represent the orientation of the tool sensors in the 
principle position. They may be calculated from Equation (C-1) assuming Frb=0: 



'jr 






0 








0^ 








0 


1 


0 


-sin(7i„) 




0 










0 


COS(7;^); 


I 0 


0 


1; 





(C-8) 

Project vector F into the gravity system using Equations (C-7) and (C-8): 



= [sin(F^)cos(0) -sin(F^)sin(<I>) cos(F^)] 



y 

- * , 



(C-9) 



6. Calculate the true formation dip and azimuth F „ from eq. (C-9) 
tan(F„) = -^ = 

sin(F^)[cos(<l>)cos(r^)siii(7;„)-sin(<l>)cos(r^)3 + co8(F^)siri(r^^)sin(r^) 
sin(F^)[cos(<I>)cos(2'^)cos(2'^) + sm(<D)sin(r^)] + cos(F„,)sin(r^)cos(r^, ) 



10 



(C-10) 



and 



^siFap ) = - sin(F„; ) cos(<D) sin(r^ ) + cos(F^, ) cos(2;^ ) 



(C-ll) 
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CLAIMS 

WHAT IS CLAIMED IS: 

1 . A method of logging of subsurface formations including a plurality of layers 
each having a horizontal resistivity and a vertical resistivity, the method 
comprising: 

(a) conveying an electromagnetic logging tool into a borehole in the 
subsurface formations; 

(b) obtaining multi-component measurements indicative of said horizontal 
and vertical resistivities using said logging tool; 

(c) defining a model of said plurality of layers, said model including a 
horizontal resistivity and a vertical resistivity associated with each of 
said plurality of layers, an azimuth angle and a inclination angle of 
said layers referenced to an earth coordinate system; 

(d) determining expected responses of the logging tool to said model; 

(e) defining a data objective fimction related to a difference between said 
e;q>ected responses and said multi-component measurements; and 

(0 using an iterative procedure including stabilization for updating the 
model to reduce said data objective fimction. 

2. The method of claim 1 wherein the stabilization comprises: 

(a) defining a model objective fimction related to changes of component 
values of said model in successive iterations; 

(b) defining a global objective function as a weighted sum of the data 
objective fimction and the model objective fimction; and 

(c) updating the model to using a sensitivity of the global objective 
fimction to model parameters. 

3. The method of claim 1 wherein said electromagnetic tool has an axis mclmed 
to a normal to said layers. 
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1 4, The method of claim 1 wherein obtaining said muhi-component 

2 measurements further comprises using at least one transmitter coil having an 

3 axis inclined to an axis of the logging tool. 

1 5. The method of claim 1 wherein obtaining the multi-component measurements 

2 further comprises using at least one receiver coil having an axis inclined to an 

3 axis of the logging tool. 

1 6. The method of claim 1 whereui obtaining the multi-component measurements 

2 further comprises using at least two receiver coils having an axis inclined to 

3 an axis of the logging tool. 

1 7. The method of claim 1 wherein said multi-component measurements comprise 

2 Hyy, and measurements. 

1 8. The method of claim 7 \^4ierein the multicomponent measurements further 

2 comprise and H^^ measurements. 

1 9. The method of claim 2 wherein the model objective function further 

2 comprises a weighted sum of individual objective functions related to said 

.3 horizontal and vertical resistivities of the plurality of layers, the dip angle and 

4 the azimuth angle. 

1 10. The method of claim 1 wherein the model objective function is related to a 

2 ratio of value of at least one component of the model after an iteration to a 

3 value of the at least one component before said iteration. 

1 11. The method of claim 1 wherein determining said expected responses further 

2 comprises: 

3 (i) obtaining an azimuth and inclination ofsaid borehole at a depth; 

4 (ii) obtaming a toolface orientation of said logging tool; 
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(iii) determining fiom said obtained borehole inclination and depth and 

said obtained toolface orientation, a relative azimuth and inclination of 
the tool to the formations. 

A method of logging of subsurface formations including a plurality of layers 
each having a horizontal resistivity and a vertical resistivity, the method 
comprising: 

(a) conveying an electromagnetic logging tool into a borehole in the 
subsurface formations; 

(b) obtaining multi-component measurements indicative of said horizontal 
and vertical resistivities using said logging tool; 

(c) defining a model of said plurality of layers, said model including a 
horizontal resistivity and a vertical resistivity associated with each of 
said plurality of layers, and a inclination angle of said layers 
referenced to an earth coordinate system; 

(d) determining expected responses of the logging tool to said model using 
a relative azimuth angle between the tool and the formation; 

(e) defining a data objective function related to a difference between said 
expected responses and said multi-component measurements; and 

(f) using an iterative procedure including stabilization for updating the 
model to reduce said data objective function. 

The method of claim 12 fiuther comprising summing and //^ components 
of said measurements and using said simuned component in the data objective 
function. 

The method of claim 12 wherein said electromagnetic tool has an axis inclined 
to a normal to said layers. 
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1 15. The method of claim 1 2 wherein obtaining said multi-component 

2 measurements further comprises using at least one transmitter coil having an 

3 axis inclined to an axis of the logging tool. 

1 1 6. The method of claim 1 2 wherein obtaining the multi-component 

2 measurements further comprises using at least one receiver coil having an axis 

3 inclined to an axis of the logging tool. 

1 1 7. A method of logging of subsurface formations including a plurality of layers 

2 each having a horizontal resistivity and a vertical resistivity, the method 

3 comprising: 

4 (a) conveying an electromagnetic loggmg tool into a borehole in the 

5 subsurface formations; 

6 (b) obtaining multi-component measurements indicative of said horizontal 

7 and vertical resistivities using said logging tool; 

8 (c) defining a model of said plurality of layers, said model including a 

9 horizontal resistivity and a vertical resistivity associated with each of 

1 0 said plurality of layers, an azimuth angle and a inclination angle of 

1 1 said layers referenced to an earth coordinate system; and 

12 (d) defining a plurality of azimuth angles and, for each of said plurality of 

13 azimuth angles, 

14 (i) inverting said multi-component measurements to give a model 

1 5 and determining a data misfit for said model, 

1 6 (ii) determining one of said plurality of azimuth angles for which 

1 7 said data misfit attains a minimum value. 
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